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Hydro-dynamic equations for a binary mixture of inelastic Maxwell models described by the Boltz- 
, ^ ' mann equation are derived. The Navier-Stokes transport coefficients are exactly obtained by solving 

, the Boltzmann equation from the Chapman-Enskog method for states close to the (local) homoge- 

neous cooling state (HCS). The knowledge of the transport coefficients allows one to analyze two 
different problems. First, we solve the linearized hydrodynamic equations around the homogeneous 
(cooling) state and identify the conditions for stability as functions of the wave vector, the dissi- 
i C| pation, and the parameters of the mixture. As happens for monocomponent systems, the analysis 

O ' shows that the HCS is unstable to long enough wavelength perturbation. As a second problem, we 

explore the validity of Onsager's reciprocal relations of a granular binary mixture. As expected, 
since a granular system is not time reversal invariant, Onsager's reciprocal relations do not apply 
for inelastic collisions. The results show that the absence of the Gibbs state (non-Maxwellian be- 
havior of the velocity distribution functions describing the HCS), the collisional cooling, and the 
occurrence of different kinetic temperatures for both species (breakdown of energy equipartition) 
are responsible for a violation of Onsager's relations. 
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I. INTRODUCTION 
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Granular fluids are usually modeled as systems composed by smooth hard spheres with inelastic collisions. The 
collisions are specified in terms of the change in relative velocity at contact but with a decrease of the magnitude of the 
normal component measured by a (positive) coefficient of restitution smaller than or equal to 1. For this interaction 
model, the Boltzmann equation (BE) has been conveniently modified to account for the inelasticity of binary collisions 
and the transport coefficients have been determined for a single gas pj as well as for multicomponent systems 0. 
However, the evaluation of those coefficients for inelastic hard spheres (IHS) is quite involved since they are given in 
terms of the solutions of linearized integral equations which can be approximately solved by considering the leading 
order in a Sonine polynomial expansion of the velocity distribution function. For this reason, other interaction models 
that simplify the mathematical structure of the Boltzmann collision integrals for IHS have been considered in the 
past few years. As for ordinary gases, the BE for inelastic Maxwell models (IMM) has been also introduced The 
IMM share with elastic Maxwell molecules the property that the collision rate is velocity independent but, on the 
other hand, their collision rules are the same as for IHS. Although these IMM do not correspond to any microscopic 
interaction potential, the cost of sacrificing physical realism is in part compensated by the amount of exact analytical 
I ■ results that can be derived from this simple model. 

Most of the studies performed for IMM have been devoted to homogeneous states, especially in the analysis of 
high-energy asymptotics of the velocity distributions. On the other hand, much less is known on the dependence of 
the transport coefficients on dissipation, especially for mixtures. In the case of a binary mixture subjected to simple 
shear flow, the relevant rheological properties of the system (shear and normal stresses) have been exactly obtained Q . 
These results show an excellent agreement with those obtained analytically for IHS in the first Sonine approximation 
and by means of Monte Carlo simulations p| . More recently, the BE for inelastic Maxwell mixtures has been solved 
from the Chapman-Enskog method for states close to the (local) homogeneous cooling state. Explicit expressions 
of the Navier-Stokes transport coefficients of IMM in d dimensions have been derived in terms of the coefficients of 
restitution and the ratios of mass, concentration, and particle sizes. Comparison with known results for IHS shows 
a reasonably good agreement, especially for the transport coefficients associated with the mass flux. 

The knowledge of the transport coefficients of IMM for mixtures allows quantitative application of the nonlinear 
hydrodynamic equations to a number of interesting problems. Here, we consider two different applications. First, 
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we study small perturbations of a spatially homogeneous state and determine the dispersion relations for the hydro- 
dynamic modes. This analysis allow us to identify the conditions for stability as functions of the wave vector, the 
dissipation, and the parameters of the mixture. As in the monocomponent case Q, linear stability analysis shows two 
transversal (shear) modes and a longitudinal (heat) mode to be unstable with respect to long wavelength excitations. 
As a second application, we explore the validity of Onsager's reciprocal relations among the different transport 
coefficients associated with the mass and heat fluxes. For ordinary gases these relations are a consequence of time 
reversal invariance of the equations of motion of the individual particles. Since a granular fluid is inherently time 
irreversible (there is a irreversible loss of kinetic energy in collisions), violation of Onsager's relations is expected 
with increasing dissipation. The interesting point here is to assess the influence of inelasticity on the failure of these 
reciprocal relations as well as to identify the origin of such a violation. 

The plan of the paper is as follows. In Sectional we give a brief summary of the Boltzmann equation for IMM and 
their corresponding balance hydrodynamic equations. The form of the Navier-Stokes hydrodynamic equations for the 
mixture is given in Sec. IIIII Sections IIVI and IVl contain to the main results of this paper. In Sec. IIVI we perform a 
stability analysis of the linearized hydrodynamic equations while Section^deals with the study of Onsager's relations 
in granular systems. We close the paper with some concluding remarks in Sec. I VII 

II. INELASTIC MAXWELL MODELS FOR A GRANULAR MIXTURE 

Let us consider a mixture of inelastic Maxwell gases at low density 0] . At this level of description, all the relevant 
information on the system is given through the velocity distribution functions / r (r, v; t) (r = 1,2, . . .) of each species, 
which obey the following set of nonlinear Boltzmann kinetic equations |4| 

(ft + v.V)/r(r,v;t)=J3 J "M/r(*). /.(*)] > (1) 

s 

where the Boltzmann collision operator J rs [v|/ r ,/ s ] is 

J r8 [vi\f r ,fs] = , ,, n / dv 2 / da 

n s {r,t)Q d J J 

x [a- 1 / r (r,v' 1 ,t)/ s (r,v^t)-/ r (r,v 1) t)/ s (r,v 2 ,t)] . (2) 

Here, m r is the mass of a particle of species r, n r is the number density of species r, u rs ^ ui sr is an effective collision 
frequency (to be chosen later) for collisions of type r-s, fid = 2-ir d ^ 2 /T(d/2) is the total solid angle in d dimensions, 
and a rs = a sr < 1 refers to the constant coefficient of restitution for collisions between particles of species r and s. 
In addition, the primes on the velocities denote the initial values {v'^v^} that lead to {vi,V2} following a binary 
collision: 

v'l = Vi - p sr (l + a" 1 ) (<T ■ gl 2 )<7, V 2 = V 2 + firs (l + O^ 1 ) (<7 • gl 2 )£ , (3) 

where gi 2 = vi — v 2 is the relative velocity of the colliding pair, <x is a unit vector directed along the centers of 
the two colliding spheres, and fj, rs = m r /{m r + m s ). The collision frequencies u> rs can be seen as free parameters in 
the model. Their dependence on the coefficients of restitution a rs can be chosen to optimize the agreement with the 
results obtained from the Boltzmann equation for IHS. 

The relevant hydrodynamic fields in a mixture are the number densities n r , the flow velocity u, and the granular 
temperature T. They are defined in terms of the distributions f r as 

dv/r(v), (4) 

pU = 2J PrUr =^2 dvm r vf r (v), (5) 
r r J 

nT = p = n r T r = £ ^ J dvV 2 f r (v), (6) 



where p r = m r n r is the mass density of species r, n = J2r n r ^ s the total number density, p = Pr ^ s the total 
mass density, V = v — u is the peculiar velocity, and p is the hydrostatic pressure. Furthermore, the third equality of 
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Eq. © defines the kinetic temperatures T r of each species, which measure their mean kinetic energies. For inelastic 
systems, in general T r ^ T so that the energy equipartition theorem does not apply. 

The collision operators conserve the particle number of each species and the total momentum, but the total energy 
is not conserved: 

j dvJrsWrJs] =0, (7) 

X] / dvm r vJ rs {v\f r , f s ] = 0, (8) 



^2 Jdv^m r V 2 J rs [v\f r J s } = ~^nT(. 



(9) 



Here £ is identified as the "cooling rate" due to inelastic collisions among all species. At a kinetic level, it is also 
convenient to introduce the "cooling rates" £ r for the partial temperatures T r . They are defined as 

Cr = = - ITT i dvm r V2j rsMfr, fs], (10) 

where the second equality defines the quantities £ rs . The total cooling rate £ can be written in terms of the partial 
cooling rates Cr a $ 

C = X>r7rCr, (H) 
r 

where x r = n r /n is the mole fraction of species r and j r = T r /T . 

The macroscopic balance equations for the mixture follow from the conditions 0-©- They are given by 

D t n r + n r V • u H — = , (12) 



Au + p- x V-P = 0, (13) 

D t T - — ^ + A (v • q + P : Vu) = -(T . (14) 

n z — ' m r an 

r 

In the above equations, D t = d t + u • V is the material derivative, 

jr = m r [dvVf r {v) (15) 
is the mass flux for species r relative to the local flow, 

P = J2 fdvm r VVf r {v) (16) 

r J 

is the total pressure tensor, and 

q = E / dvim r ^ 2 V/ r (v) (17) 

r •* 

is the total heat flux. The balance equations i|12|) - (|14[l apply regardless of the details of the model considered for 
inelastic collisions. However, the influence of the collision model appears through the dependence of the cooling rate 
and the hydrodynamic fluxes on the coefficients of restitution. 
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As happens for elastic collisions the main advantage of using IMM is that a velocity moment of order k of the 
Boltzmann collision operator only involves moments of order less than or equal to k |9j . This allows one to determine 
the Boltzmann collisional moments without the explicit knowledge of the velocity distribution functions. The first 
few moments of the Boltzmann collision operator J rs [fr,fs] have been explicitly evaluated in Ref. 0. In particular, 
the cooling rates Crs are given by 



_ 2uj rs 

Srs j psr\^ t Oi rs ) 



Psr M , ,0 r + S fl sr (l + a rs ) - 1 . 

2 S dp s p r 



(18) 



where p r — n r T r and 



= ^5>r*. d9) 



Equation (|18|) can be used to fix the parameters Lu rs . The most natural choice to optimize the agreement with the 
IHS results is to adjust the cooling rates £ rs for IMM, Eq. (|18fl . to be the same as the ones found for IHS. Given 
that the coolin g ra tes are not exactly known for IHS, one can estimate them by considering their local equilibrium 



Imgra 

approximation |10| . With this choice, the collision frequencies u> rs are given by 



'TV 



" '<-< (-^) -o, (20) 



where o~ rs — (a r + cr s )/2, a r being the diameter of particles of species r. In addition, vq — (2T ^ r m" 1 ) 1 / 2 1S a 
thermal velocity defined in terms of the global temperature T. Upon deriving i|2(J|l use has been made of the fact that 
the mass flux j r vanishes in the local equilibrium approximation. In the remainder of this paper, we will take the 
choice (|2"n)) for uj rs . 

III. NAVIER-STOKES HYDRODYNAMIC EQUATIONS 

Needless to say, the usefulness of the balance equations (|12ll - H14|) is limited without further specification of the 
fluxes and the cooling rate on space and time. However, for sufficiently large space and time scales, one expects that 
the system achieves a hydrodynamic regime in which all the space and time dependence of the distribution function 
occurs through a functional dependence on the hydrodynamic fields. This functional dependence is made local in 
space and time by writing / r (v) as a series expansion in powers of the gradients of the hydrodynamic fields. This 
special solution is called a normal solution and can be obtained by applying the Chapman- Enskog method |Tlj to the 
BE. One important difference with respect to the conventional Chapman-Enskog method for ordinary gases is that 
the reference state (zeroth-order approximation) of the expansion is not the local equilibrium state but the so-called 
(local) homogeneous cooling state, whose explicit form is not known ^(j ■ Very recently, the Chapman-Enskog solution 
to the BE has been worked out to first order in gradients of the fields. The corresponding constitutive equations 
found up to this order for the binary mixture are [(| 

= _™^ Vxi _ P DpVp _ LMT t h = -jr, (21) 
P P T 

q = -T 2 D"V Xl - nVp ~ AVT, (22) 



Pm = pSke - V (\7eu k + VfeMf - • u^j . (23) 

The transport coefficients are the diffusion coefficient D, the thermal diffusion coefficient D', the pressure diffusion 
coefficient Dp, the Dufour coefficient D" , the thermal conductivity A, the pressure energy coefficient k and the shear 
viscosity coefficient 77. In contrast to the results previously derived for IHS 0, all the above transport coefficients 
have been exactly obtained in terms of the coefficients of restitution and the parameters of the mixture (masses, sizes, 
and composition). Technical details of the Chapman-Enskog solution to the BE (JIJ can be found in Ref. ,6j. Here, 
for the sake of completeness, the explicit expressions of the Navier-Stokes transport coefficients as well as the cooling 
rate C are displayed in the Appendix . 
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The expressions for the mass flux (|21() , the heat flux ill' 21) , the pressure tensor (|23ll , and the cooling rate (|A6|) provide 
the necessary constitutive equations to convert the balance equations (|T^|l - ifTH) into a closed set of six independent 
equations for the hydrodynamic fields. Since the irreversible fluxes have been represented in terms of the gradients of 
the mole fraction x\, the pressure p, the temperature T, and the flow velocity u, it is convenient to use these as the 
independent hydrodynamic variables. This means that, apart from the balance equations (|13|l and l|14|) for u and T , 
respectively, we also need the corresponding balance equations for x\ and p. These equations can be easily obtained 
from I112[I and l|14|l and are given by 

D t x 1 + - r ^ V-ji=0, (24) 

2 

D t p+pV-u+ — (V-q+P:Vu) = -Cp. (25) 
net 

Therefore, when the expressions of the fluxes and the cooling rate £ are substituted into the balance equations JT3J, 
(|14p. I|24l) . and l|25|l one gets a closed set of hydrodynamic equations for x±, u, T, and p. These are the Navicr-Stokes 
equations for the granular binary mixture. They are given by 

D tXl = -^—^ • ( V -^DV X1 + P -D p Vp + Id'Vt) , (26) 
rmim2 \ P p i J 

(D t + C)P+ —J-P^ ' u = \ V ■ ( t2d "Vxi + «Vp + AVT) 

+ (\7 e u k + V k ut - ~5|«V • u^j Viu k , (27) 

(D t + ()T+ -j-pV ■ u = -V • DVxi + -D p Vp + ^D'VT 

an n mim 2 VP P -L J 

+-I-V ■ (T 2 D"Vxi + k\7 p + AVT) 
an 

+ J^i? ^V^iifc + Vfcii* - ~£a*V -u^V^Ufc, (28) 

+ P~ Y ^iP = P _1 V fe 77 ^V^Ufc + V fe U£ - ^(5«V ■ u^j . (29) 

IV. LINEAR STABILITY ANALYSIS OF THE HYDRODYNAMIC EQUATIONS 

As said in the Introduction, one of the simplest application of the hydrodynamic equations is a stability analysis of 
the nonlinear hydrodynamic equations (|26l) - (|29|l with respect to the homogeneous state for small initial perturbations. 
The linearization about the homogeneous solution yields partial differential equations with coefficients that are inde- 
pendent of space but depend on time since the reference (homogeneous) state is cooling. As in the monocomponent 
case pj, this time dependence can be eliminated through a change in the time and space variables, and a scaling of 
the hydrodynamic fields. 

Let us assume that the deviations dy a (r,t) — y a (r,t) — yna{t) are small. Here, 6y a (r,t) denotes the deviation of 
{x\,u, T, p} from their values in the homogeneous state, which is indicated by the subscript H. We introduce the 
following dimcnsionless space and time variables: 

r= fdt'v B {t'), r' = ^|r, (30) 
Jo v 0H (t) 



where vh{^) = {^d/^V : ^) n HCi2 v 0H is & n effective collision frequency and vqh — •\/2Tjj(mi -\- m^) J m\m2- Since 
{xih,uh,Th,Ph} are evaluated in the homogeneous cooling state, then 

d t xiH = 0, u fl = 0, d t mTff = d t \np H = -( H . (31) 



A set of Fourier transformed dimensionless variables are then introduced as 



5xik{T) 5u k (r) 
Pk{r) = , w k (r) - 



XlH Vqh(t) 

where Sy^a = {Sxi^, <5iik, <5Tk, <5pk} is defined as 



<% ka (k,r) = /" dr' e - lk r '^ Q (r', r). 



(32) 



(33) 



Note that here the wave vector k is dimensionless. In terms of the above variables, the transverse velocity components 
Wki = Wk— (wk-k)k (orthogonal to the wave vector k) decouple from the other four modes and hence can be obtained 
more easily. They obey the equation 



where = (h /vh is given by Eq. (|A6|I and 



v H r] 



PhvI h ' 



The solution for Wki (k, r) reads 



where 



w k j_(k,r) = w kJ _(0) exp[s j. (/c)r], 



sx(k) = ~C-V*k 2 . 



(34) 

(35) 

(36) 
(37) 



This identifies two shear (transversal) modes. We see from Eq. (|37() that there exists a critical wave number k^_ given 
by 



2rf 



1/2 



(38) 



This critical value separates two regimes: shear modes with k > fc^ always decay while those with k < k c ^ grow 
exponentially. 

The remaining modes are called longitudinal modes. They correspond to the set {pk,#k,IIk} along with the 
longitudinal velocity component u>k|| = Wk • k (parallel to k). These modes are the solutions of the linear equation 



8t 



(MW+ikMi$ + tfM$>)6yw(T), 



(39) 



where 5yk a (r) denotes now the four variables {pk, #k, Ilk, K>k|| }• The matrices in Eq. I|39|l are given by 



/ 















-x Y {dC/dxi) Tp 


1 ^* 
2^ 


-C 


°) 




-x Y (dC/dxi) TiP 


1 A* 
2 C » 


-c 





V 











\C 1 



(40) 



M« = 



/0 
' 








y o o i<Vi±±<4 



_ 2 

d 
d+2 

d 





(41) 
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In these equations, /i = rai/m2, and we have introduced the reduced Navier-Stokes transport coefficients 



D* 



D* 

2 ' W 



mim 2 n 2 H v 2 H ' 



m 1 m 2 n 2 H vl H ' 



(43) 



„ = uhTrD" 
n H v 2 H 



U 0H 



A* = 



(44) 



The longitudinal modes have the form exp[s„(fc)r] for n = 1, 2,3,4, where s n (fc) are the eigenvalues of the matrix 
] = M^ ) + ikM^ 1 ' + k 2 M^ 2 \ namely, they are the solutions of the quartic equation 



detlM-slI =0. 



(45) 



The solution to l|45|) for arbitrary values of k is quite complex. Interestingly, the expression of M through first order in 
k does not depend on the assumed forms of the mass flux, the pressure tensor, and the heat flux. As a consequence, 
the corresponding eigenvalues can be considered as known exactly up through this order. They are found to be the 
eigenvalues of the matrix of M (°) : 



(0,0, 



(46) 



Hence, at asymptotically long wavelengths (k = 0) the spectrum of the linearized hydrodynamic equations is comprised 
of a decaying mode at — £*/2, a two-fold degenerate mode at 0, and a three- fold degenerate unstable mode at (* /2. 
For finite k, the longitudinal modes can be calculated for small k by a perturbation expansion: 



Sn (fc) = s (o) + fc s W + fe^) + ... 



where is given by Eq. (|46() . Through order k 2 , the coefficients (for n = 1, 2, 3, 4) are 

s«=0, n= 1,2,3,4, 



(47) 



(48) 



,(2) 



,(2) 



1 1 



dlnC 



(D' : 



2 X\fJL + x 2 



D* 
1 



2C* (xi/i + ar 2 )(l+/x) 



(2jD '* +jD *)_ D 



(49) 



,(2) 



ainc 

9a; i 



d+1 



1 — n 

XlfX + X 2 



(D'*+D* p ) 



d(* {xin + x 2 )(l+n) rf (A * + K * )j 



(50) 



s ( 2 ) = -2^i«* - - M f51) 

4 d ' rfC* (a;i/i + a;a)(l+M)' 

Since the Navier-Stokes order only applies through order k 2 , the solutions l|46|l and H48fl -151 [) are relevant to the same 
order. 

To illustrate the general dependence of the hydrodynamic modes on the reduced wavenumber k, in Fig. ^ we show 
the dispersion relations for the case an = a 22 = ol\ 2 = a = 0.9, <J\ja 2 = 1, x\ = 0.2, and m\/m 2 = 4. From top 
to bottom the curves correspond to the shear modes, given by Q37f) . and the remaining four longitudinal modes. We 
only represent the real parts of the hydrodynamic modes. Furthermore, we have also included the results obtained 
in the case of IHS by using the leading Sonine approximation for the transport coefficients It is seen that the 
shear mode and one of the longitudinal modes ("heat" mode) are positive for k < k\ and k < respectively. Here, 
k^ is given by Eq. I|38|l while the critical value fc£ for the heat mode may be found from the condition dct |M| = 0. 
Therefore, initial long wavelength perturbations of the homogeneous cooling state that excite shear and heat modes 
will grow exponentially, representing an instability of the reference state. In addition, we also observe that the results 
for IMM compare quite well with those given for IHS showing again the reliability of IMM as a toy model to capture 
the main trends observed for granular fluids. 
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FIG. 1: Dispersion relations in the three-dimensional case for an = «22 = ai2 = a = 0.9, ai = oi, xi = 0.2, and mi/m.2 = 4. 
From top to bottom the curves correspond to the shear modes and the remaining four longitudinal modes. The solid lines are 
the results derived for IMM while the dashed lines refer to the results for IHS. 



ONSAGER'S RECIPROCAL RELATIONS 



As a second application of the explicit knowledge of the Navier-Stokes transport coefficients, in this Section we 
explore the validity of Onsager's reciprocal relations for inelastic Maxwell mixtures. In the case of elastic collisions, 
Onsager's relations establish symmetry properties between the Onsager phenomenological coefficients L sr , L sq , 
L qq , and L qs associated with the mass and heat fluxes of a gas mixture. These coefficients are defined through the 



(linear) constitutive equations |JJ 



J>2 ' 



(52) 



d + 2 



if. 



-L qq VT 



2 k 



V T 



Here, we have introduced the gradient of chemical potential per unit mass (p, s ) given by 



Onsager's relations state Q 



T 



Vln(x s p). 



(53) 



(54) 



(55) 



Relations (|55|l are a consequence of time reversal invariance of the equations of motion of the individual particles. 
Since the main feature of a granular gas is the irreversible loss of kinetic energy in collisions, one expects relations 155|) 
not to apply when a rs ^ 1. Now we want to assess the effect of dissipation on the violation of Onsager's theorem. 

To identify the Onsager phenomenological coefficients from the constitutive equations l|52|) and i|53|) . Eqs. (|21|) and 
(12211 must be rewritten in terms of the gradients of chemical potential p, r and temperature T . To do that, first note 
that Eq. I|54|) leads to the relation 



P1P2 (V/^i)t - (V/x 2 )t «4«-2 



np 



T 



np 



(ni2 — mi)V hip, 



(56) 



where use has been made of the identity Vx\ = —Vx2- When Eq. (|56|l is substituted into Eqs. I|21J) and (|22(l . one 
gets the following expressions for the mass flux ji and the heat flux 3 q defined in Eq. H53JI : 



Ji 



m 1 m 2 piP2 >i)t 



(VAt 2 ) 5 



C p Vlnp-^'VT, 



(57) 
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FIG. 2: Plot of the reduced function P(a) as a function of the coefficient of restitution a in the three-dimensional case for 
cri = (72, Xi = 0.2, and three different values of the mass ratio: (a) mi/m2 = 0.5, (b) m\/rri2 = 2, and (c) mi/7712 = 4. 



P1P2D 



( P1P2T 2 „ _ d + 2 T m 2 - mi 
" ^ np 2 p 2 

(, d + 2 m 2 — mi ,\ 

-C p V lnp - A - — ? l -pD' VT, 

F \ 2 m\m 2 



(V^i)t - (V^ 2 )r 



T 



(58) 



where 



and 



C p = p£> p - ^p-{m 2 - mi)D, 



C p =pn 



d+2 rTl m 2 -mi„ n x n 2 2 

^ P ^— T (m 2 -mi)D . 



mim 2 



(59) 



(60) 



The Onsager phenomenological coefficients can be easily identified from Eqs. i|57[) and (|58|l when one takes into 
account the constitutive equations (|52|l and l|53|l . They are given by 



mim 2 pip 2 

Mi = -^12 = 5 i*) = plJ-> , 



(61) 



-I 



92 



T 2 p lP2 
np 



D 



d+2 Tpip 2 
~2 



(m 2 -mi)D, L qq = A 



<i + 2 m 2 - mi , 



zr. 



mim 2 



(62) 



The first Onsager relation L\ 2 — L 2 \ is trivially verified since, according to Eq. I|A3(1 . the diffusion coefficient is 
symmetric under the exchange 1 <-> 2. The second Onsager relation requires that L\ q — L q ±. In addition, according to 
Eqs. l|57|) and there are two new contributions proportional to Vp not present in the structure given by Eqs. l(S^|) 
and (|53[) for the fluxes ji and 3 q , respectively. Consequently, Onsager's relations would require that these contributions 
should vanish for any value of a rs , i.e., C p — C' p = 0. In terms of the reduced coefficients {D* , D*, D'* , D"* , k* , X*} 
defined in the Appendix, the conditions L± q = L q \, C p = 0, and C' p — lead, respectively, to the following conditions 



P({a rs })=0, Q({a rs }) = 0, R{{a rs }) = 0, 



(63) 



where 



P({a rs }) = D" 



d + 2 1 - p 2 
~~2 ~p~ 



D* 



1 + p x 2 + px x 

P X\X 2 



(64) 



Q({a rs }) =Dt- x x x 2 - D* , 
x 2 + pxi 



(65) 
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R({a rs }) = K* - d+ n 21 tl2 Q-x lX2 \ » D"*. (66) 
2 fi x 2 + [JiXi 

In summary, if Onsager's relations held, then the reduced functions P, Q, and R would vanish for any value of the 
parameters of the mixture. In the elastic limit (a rs = 1), our results actually show that -P(l) = <3(1) = R{1) = 0, 
so that, Onsager's relations are exactly verified. Nevertheless, for inelastic collisions, the functions P, Q, and R are 
different from zero, as expected. There are mainly three independent sources of discrepancy: (i) non-Gaussianity of 
the distribution functions of the homogeneous cooling state, which is measured through the cumulants c r and are 
given by Eq. (|A20|) ; (ii) energy nonequipartition as measured by the deviation of 7 from unity, and (iii) time evolution 
of the granular temperature, which is accounted for by the cooling rate £. To illustrate the influence of dissipation 
on the breakdown of Onsager's symmetry relations, the functions P, Q, and R are plotted in Figs. [3J |3 and 
respectively, as functions of the coefficient of restitution a. Here, for the sake of simplicity, we have assumed again a 
common coefficient of restitution (a rs = a) for d — 3 with X\ — 0.2, o\ — 02, and different values of the mass ratio 
7711/7712. As might be expected, the deviation of P, Q, and R from zero increases with decreasing a, and in general is 
more significant as the mass disparity increases. In addition, while the magnitude of Q(a) is quite small, the function 
i?(a) grows very fast with dissipation. The validity of the pure Onsager's relation L\ q = L q \ is tested by the function 
P{cx). Except for very weak dissipation, the violation of this reciprocal relation is quite important, especially in the 
case toi/to2 = 4. 
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VI. CONCLUDING REMARKS 

In this paper we have derived the hydrodynamic equations for a granular binary mixture from the Boltzmann 
kinetic theory for inelastic Maxwell models (IMM). In the Boltzmann equation for IMM the collision rate of inelastic 
hard spheres (IHS) is replaced by an effective collision rate independent of the relative velocity of the two colliding 
particles. This simplification allows one to compute the velocity moments of the Boltzmann collision integrals without 
the explicit knowledge of the distribution functions. Thanks to this property, the Navier-Stokes transport coefficients 
of the mixture and the cooling rate have been exactly obtained in terms of dissipation and the parameters of the 
mixture (masses, sizes, composition) [(|. These results contrast with the ones previously derived for IHS, where the 
transport coefficients were approximately determined by using the leading terms in a Sonine polynomial expansion of 
the distribution functions. 

The explicit knowledge of the transport coefficients and the cooling rate allows one to study two different prob- 
lems. As a first application, we have performed an analysis of the linearized hydrodynamic equations around the 
homogeneous cooling state. As in the monocomponent case for IHS [lj, our stability analysis shows that the homo- 
geneous cooling state is unstable to long enough wavelength perturbations and consequently becomes inhomogencous 
for long times. Specifically, there are two shear modes and a longitudinal mode which become unstable for small 
values of the wavenumber. Therefore, small perturbations or fluctuations around the homogeneous state that excite 
the above modes will grow exponentially As a second application of a hydrodynamic description, we have explored 
the breakdown of the Onsager relations between transport coefficients associated with the mass and heat fluxes of 
a granular mixture. For elastic systems, Onsager's reciprocity relations are a consequence of the fact that the me- 
chanical equations of motion (classical as well as quantum mechanical) are symmetric with respect to time inversion. 
This microscopic property leads to Onsager's theorem. Given that the above time reversal invariance is broken in 
dissipative systems, violation of Onsager's theorem is expected for inelastic collisions. Here, we have analyzed the 
effect of the dissipation on such a violation. Our analysis shows that violation of Onsager's relations in granular gases 
has basically three distinct origins. First, the deviation of the homogeneous cooling state from the Gibbs state is 
responsible for the coefficients C\ and C2 [defined by Eq. (|A20|l ] being different from zero. At a quantitative level, 
this effect is relatively quite small. Second, the effect of collisional cooling occurs through the presence of the cooling 
rate £• Finally, the effect of different partial temperatures is expressed by the factors 7 r and 8 r in the different terms 
involved in the calculation of the transport coefficients. Each one of these effects is a different reflection of dissipation 
in collisions. 

In summary, there is growing theoretical support for the usefulness of a hydrodynamic description for granular 
systems under rapid flow conditions. However, in spite of the similarities between granular and normal fluids, the 
extension of properties of ordinary fluids to those with inelastic collisions must be carried out with caution. As shown 
here, the homogeneous (cooling) state is unstable to long wavelength perturbations and the familiar Onsager relations 
do not apply. These are some examples which makes these systems quite different from gases of elastic particles, like 
molecular gases. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF THE TRANSPORT COEFFICIENTS 

In this Appendix we display the expressions of the transport coefficients {D, D p , D' , D" , k, A, 77} defining the mass 
flux (|21|l . the heat flux (|22|l . and the pressure tensor (|23|> . Let us consider first the coefficients associated with the 
mass flux. The expressions of the coefficients D', D pi and D are given, respectively, by Q 



(Al) 
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D = 



PT 



9X X , T 



dC 

9xi / , 



1 p Wp 



" ~2 C 



where 



d-l 

= ^=ncr 12 v Q , v 



2T 



mi + to 2 



(A3) 



(A4) 



4-^/7? " V 77117712 

In these equations, /i = mi/f)i2 is the mass ratio, £* = C/^Oj and we have introduced the (reduced) collision frequency 



- d (x 2 ^i+xi^ 12 ) I — 



1/2 



(1 + aia), 



(A5) 



where d\ = l/(/i2i7i) and 6 2 = l/(/Ji272)- The temperature ratio 7 = T\/T 2 is determined from the condition 
£1 = C2 = C I n reduced units, the cooling rate = C,i/ v o is given by 



2^2 /c7i\ d 1 1/2 2 8 



+ 2 



'1P2 



1/2 



(l + ai 2 ) 



M21 , 
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l-^(l + ai 2 )- 



+ 6*2 



(A6) 



while the expression of Q = C$/vq can be easily obtained from l|A6() by the exchange 1 <-> 2. It must be noted that £* 
depends explicitly on xi and also implicitly through its dependence on the temperature ratio 7. Since ji = — j 2 and 
Vxi = — Va; 2 , it is expected that D is symmetric with respect to exchange of particles 1 and 2 while D p and D' are 
antisymmetric. This can be easily verified by noting that X171 + £272 = 1. 

The case of the heat flux is more involved. The transport coefficients D", k, and A are given by 



D" = D'l + D'2, k = Ki + k 2 , A = Ai + A 2 



where 
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nT 



(mi +m 2 )v r ' (to 1 + 7772)^0 r (toi + 7772)^0 

By using matrix notation, the coupled set of six equations for the unknowns 

{DT.DZ*. a;, A5> 

can be written as 

■^(Tfj'A^fj' . 

Here, X a i is the column matrix defined by the set 1|A9|1 and A aa i is the matrix 
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A*21 



{/z 2 i(l + a 12 ) [d + 8 - 3^2i(l + ««)] - 3(d + 2)} , 



(A7) 
(A8) 

(A9) 
(A10) 



(All) 



(A12) 



13 



,., , h+0 2 \ 1/2 2 (1 + « 12 ) 3 

v 12 = -12x1 [ a a Pi2P 2 l" 



d{d + 2) 



The elements of the column matrix Y are 
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Here, we have introduced the cumulants c r measuring the deviations of the zeroth-order distribution functions fr 
from a Maxwellian, 
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The reduced quantities A* s are given by 



A 12 = k X 
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The corresponding expressions for A^, ^22, and 1/21 can be inferred from Eqs. (|A1 2|) . (|A13jl . and (|A2 1|) by exchanging 

l<->2. 

The solution to Eq. (|A10|) is 



(A22) 



This relation provides an explicit expression for the coefficients D"* , n* , and A* in terms of the coefficients of restitution 
and the parameters of the mixture. 

Finally, the shear viscosity coefficient rj can be written as 



V= — {Vi +V2) , 
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where the expression of the (reduced) partial contributions -q* is 

^i7i( 2t 22 - (*) - 2ar 2 72Ti2 
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The quantities and r* 2 are given by 
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d(d + 2) 



jU 2 i(l + "12) / 6>i + 6*2 



nf2 



1/2 r 



_ ^21(1 + 012) 
d + 2 



(A24) 



(A25) 



8^2 ( 1 +B 2 \ 1/2 p x 2 2 
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A similar expression can be obtained for by just making the exchanges 1 «-> 2 
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